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Abstract. Deep, near-infrared imaging surveys have been motivated by the desire 
to study the rest-frame optical properties and stellar content of galaxies at high 
redshift. Here we briefly review their history, and illustrate one application, using 
HST NICMOS imaging of the Hubble Deep Field North to examine the rest-frame 
optical luminosities and colors of galaxies at 2 < z < 3, and to constrain their 
stellar masses. The rest-frame _B-band luminosity density at z ~ 2.5 is similar to 
that in the local universe, but the galaxies are evidently less massive, with rapid 
star formation and low mass-to-light ratios. There are few candidates for red, non- 
star-forming galaxies at these redshifts to the HDF/NICMOS limits. We estimate 
a stellar mass density at 2 < z < 3 that is ~ 5% of the present-day value, with an 
upper bound of < 30%. Future headway will come from wide-field ground-based 
surveys, multiplexing infrared spectrographs, and new space-based facilities such 
as HST/WFC3 and SIRTF/IRAC. 

1 Introduction and history 

Observing proposals and article introductions almost universally list a set 
of basic, interrelated themes that motivate deep near-infrared (NIR, here 
regarded as 1-3 jum) blank sky surveys: 

• The integrated stellar spectra of normal galaxies peak in the NIR 

• The stellar component of the extragalactic background peaks in the NIR 

• NIR light measures familiar rest-frame optical wavelengths at high z 

• Optical/NIR rest-frame light comes primarily from mid- to low-mass 
stars with long lifetimes relative to Hq 1 

• NIR light traces total stellar content/mass 

• Evolutionary corrections are smaller and easier to model 

• fc-corrections are small or negative and similar for most galaxy types 

• Effects of dust extinction are greatly reduced 

• Access to z > 6, where galaxy light shifts beyond optical wavelengths. 

The earliest, heroic efforts jl]j^] used single-element photometers to search for 
sky background fluctuations from primeval galaxies (PGs) whose light might 
be redshifted beyond the optical wavelength range. The field really came to 
life with the advent of array detectors, leading to the first faint NIR imaging 
surveys [gJQ]. It is interesting to read these and other early papers, where we 
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find discoveries, concerns and hypotheses that have stayed with us ever since: 
PGs, extremely red objects ("EROs" as PG candidates, or high-z ellipticals, 
or dust-enshrouded galaxies), extremely blue objects (with rapid, cosmologi- 
cally significant star formation), ERO clustering, UV-excess ellipticals, pho- 
tometric redshifts, NIR number counts (to constrain space curvature and/or 
galaxy evolution), etc. While many of the issues remain the same today, the 
data quality has advanced tremendously, largely driven by progress in array 
technology, and most recently by the leap into space with iJST/NICMOS 
(affording high angular resolution and far lower backgrounds). Survey sensi- 
tivities have improved by a factor of ~ 1000, source densities have increased 
~ 200-fold, and we can now image the detailed morphologies of high redshift 
galaxies in their rest-frame optical light. 



2 Infrared observations of the Hubble Deep Fields 

The HST WFPC2 and STIS observations of the Hubble Deep Fields (HDFs, 
North and South) are the deepest optical images of the sky, and correspond- 
ingly deep NIR observations of these areas are valuable for all the reasons 
outlined above. The HDF-N was observed from the ground in several differ- 
ent NIR programs [|||!,Q,|§||| , while the HDF-S has been imaged from the 
ESO NTT Jh| and more recently to with ISAAC on the VLT fy]. The deep- 
est observations at 1.1-1.6 /im have come from HST NICMOS imaging of 
the HDF-N |l2],[l3) and for the HDF-S NICMOS field (w) (which is distinct 
from the HDF-S WFPC2 field). The depth and angular resolution of our 
"wide-field" (only ~ 6 arcmin 2 , smaller than the first NIR array surveys!) 
HDF-N/NICMOS program, combined with the great wealth of supporting 
imaging and spectroscopy at other wavelengths from the ground and from 
space, make this a premier resource for studying the NIR properties of galax- 
ies at high redshift. Discussion of the NIR morphological and photometric 
properties of galaxies at < z < 2 JlBf and 2 < z < 3.5 [0, and of galaxy 
candidates at z ^> 4 |l6|,p7|Jl^,[l^,p0 21 have appeared elsewhere. 



3 Stellar populations of galaxies at 2 < z < 3 

We have carried out a detailed study p2| of the stellar population proper- 
ties of star-forming "Lyman break galaxies" (LBGs) from the HDF-N at 
2 < z < 3.5, using NICMOS data to extend previous work based on ground- 
based NIR photometry |^3| . Using a sample of 33 spectroscopically confirmed 
galaxies [|4|, we compared 7-band (0.3-2.2 ^m, observed frame) photome- 
try to empirical spectral templates for nearby galaxies and to population 
synthesis models in order to evaluate constraints on the LBGs' stellar con- 
tent and evolutionary histories. The LBGs are much bluer than local, Hubble 
sequence galaxies, and than comparably luminous HDF galaxies at lower red- 
shift (Fig. [I]) , but are similar to nearby, UV-bright starburst galaxies Jl^,^2| . 
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HDF "100 brightest galaxies" sample (L > L* B (z = 0)) 
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Fig. 1. Rest-frame U — B color vs. redshift for HDF-N galaxies with rest-frame 
M B < -20.3 {Q M = 0.3, Q A = 0.7, h = 0.7). Filled and open symbols indi- 
cate galaxies with spectroscopic and photometric redshifts, respectively. There is 
a strong color trend with redshift; nearly all galaxies at 2 < z < 3 are bluer than 
comparably luminous galaxies at z < 1. One redder object with Zphot = 2.7 is 
marked; this object is also a radio and X-ray source J26| . and may be detected 
at 15 /im f27[ and 1.3 mm [Eq] as well. 



Even with high-quality HST optical/infrared photometry, we find only 
weak constraints on most parameters of the LBG stellar populations, with 
degeneracies between age, star formation time scale, metallicity and extinc- 
tion. Perhaps the best constraints, however, are those on the total stellar 
mass M. (see Fig. ||). If the LBG star formation history ^(t) is modeled by 
an exponential, & oc e - */ r , with t and r (and extinction) as free parameters, 
then with fixed assumptions about the IMF and metallicity, the typical 68% 
confidence interval on log M. is approximately ±0.25 dex. For LBGs with L* 
UV luminosities (29), the inferred stellar masses (assuming a Salpeter IMF, 
and varying the model metallicities) are 1 to 2 x \Q w M.q for Qm = 0.3, 
Qa — 0.7, h = 0.7 (used here unless otherwise noted). These are ~ l/10th 
the stellar masses of L* galaxies today |50|. We may compare these masses 
to those derived from virial estimates using nebular line-widths and HST- 
measured sizes plp2] |, which are also ~ 10 lo .M©. This suggests that these 
kinematic measurements underestimate the full mass of the dark matter halo, 
reinforcing a point emphasized by Max Pettini in this volume. 
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Fig. 2. Stellar masses for HDF-N Lyman break galaxies derived from stellar pop- 
ulation model fitting |^2| . The filled points show best-fitting mass estimates using 
solar metallicity, Salpeter IMF models, and assuming exponential star formation 
histories, while the error bars show 68% confidence intervals. The downward tri- 
angles are upper mass limits allowing for the presence of an underlying, maximum 
Ai/L stellar population formed at z = oo. The dashed line shows the present-day 
characteristic M* stellar mass Ipotl for the same IMF. 



We set an upper bound on the allowable stellar mass by considering how 
much light from a hypothetical, maximally old stellar population (formed 
at z = oo) could be hidden beneath the glare of the young, star-forming 
population. On average, this upper bound is a factor of ~ 5 to 6x larger 
than the mass derived for the "young" models. If this were generally the 
case, however, then virtually all galaxies with previous generations of star 
formation at z > 3 must also be forming stars rapidly at 2 < z < 3. We see 
very few candidates for mature, red, non-star-forming galaxies in this redshift 
range, even with a NICMOS-selected sample where photometric redshifts 
should, in principle, readily identify such galaxies if they are present (see 
Fig- 0). 
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4 The optical luminosity function and stellar mass 
density at 2 < z < 3 

Using photometric and spectroscopic redshifts and NICMOS photometry for 
HDF-N galaxies, we may examine the rest-frame B-band luminosity dis- 
tribution of galaxies at 2 < z < 3 (Fig. ||). Over the luminosity range we 
can examine, this is not dissimilar to the local B-band luminosity function 
(LF) Q. Galaxies with 2 < z phot < 3 and M B < -18.75 contribute a to- 
tal blue luminosity density ps = 5.2 x 10 26 erg s _1 Mpc~ 3 . Without further 
correction for incompleteness or extrapolation to fainter magnitudes, this is 
~ 1.5 x that from the integrated 2dF LF, and nearly equal to that from the 
preliminary SDSS LF @. 




M B (AB) 

Fig. 3. Rest-frame B-band luminosity distribution for galaxies at 2 < z < 3 in 
the HDF-N. Horizontal bars indicate magnitude ranges where spectroscopic or 
photometric redshifts, or a mix of the two, have been used. The data have not 
been corrected for incompleteness, which may affect the fainter points. Error bars 
indicate Poisson uncertainties only. The local luminosity function from the 2dF 
survey |3^] is shown for comparison. 

Although the optical luminosity densities are similar at z = and z = 3, 
the implied stellar mass densities are quite different. Galaxies at 2 < z < 3 
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are far bluer than local counterparts (Fig. [I]), indicating much smaller mass- 
to-light ratios. From our modeling |2j for 22 HDF-N LBGs at 2 < z < 3, 
we derive an average (Ad /Lb) = 0.10 to 0.25 (in solar units), depending 
on assumptions about metallicity and IMF. These values refer to emergent 
luminosities, i.e., here Lb is not corrected for the effects of extinction - the 
intrinsic Ad /L for the stellar populations are smaller. There is a trend of Ad /L 
with rest-frame color, amounting to a factor of ~ 2 over the observed color 
range of the spectroscopic LBG sample, but we neglect this here and apply 
(Ad/ Lb) from the spectroscopic LBGs to the photometric redshift sample. 
The majority of galaxies with photometric redshifts 2 < z p hot < 3 have colors 
similar to or bluer than those in the spectroscopic sample (Fig. 0). 

Restricting our attention to a Salpeter IMF, the starburst dust attenu- 
ation law poj ], and model metallicities 0.2 to IxZq, we estimate a stellar 
mass density = 1.7 to 2.9 x lO 7 7W0Mpc -3 at these redshifts. Comparing 
this to the present-day stellar mass density computed from the 2dF+2MASS 
K— band luminosity function |]30| , using the same IMF assumptions, we find 
p*(z = 2.5)/p*(z = 0) = 0.034 to 0.058. It is also 8-14x smaller than the 
estimated mass density from bright galaxies at z rj 0.9 [^6| . This is presum- 
ably a lower limit for several reasons. First, our assumed cosmology results 
in a nearly minimal luminosity density for currently acceptable values of the 
cosmological parameters. An Einstein-de Sitter model increases by the lu- 
minosity and mass densities by ~80%. Second, we have made no corrections 
for incompleteness in the NICMOS-selected galaxy sample, nor any attempt 
to extrapolate to objects fainter than the HDF/NICMOS detection limit. 
Finally, as described in §2, the galaxy masses may be larger if there were 
earlier generations of star formation. If we assign every LBG its maximum 
(at 68% confidence) stellar mass (see §2) allowing for a older generation of 
stars formed at z = oo, then for the adopted cosmology we may set a con- 
servative upper bound on the total stellar mass density contained within 
NICMOS-detected galaxies, p*(z = 2.5) /p^z = 0) < 0.30. 

This upper bound is barely consistent with the hypothesis that all stars 
in present-day galactic spheroids formed at z > 3 j37|. In this scenario, all 
galaxies must have already formed most of their stars at z > 3, but must 
also be forming more stars at z k 2.5 (since there are are few candidates 
for evolved, non-star-forming HDF galaxies at that redshift). As described 
above, a more direct accounting for the mass present in NICMOS-detected 
galaxies at z w 2.5 implies a much smaller fraction, ~ 5%, of the present-day 
stellar mass density. 

5 Future directions 

Infrared surveys described here have really just begun to address the most 
important questions about the mass assembly history of galaxies. A new gen- 
eration of infrared instrumentation, on the ground and in space, will carry 
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us further along down this road. Deep NIR surveys are still limited to very 
small solid angles and thus volumes at high redshift. New large format de- 
tectors will greatly improve this situation, although unfortunately very few 
cameras for 8-10m telescopes (whose aperture is really needed to study z > 2 
galaxies) are being configured with wide fields of view. The first multi-object 
NIR spectrographs are just now coming on line; these will permit wholesale 
spectroscopy of distant galaxies at rest-frame optical wavelengths, offering a 
means of measuring kinematic masses and chemical abundances at high red- 
shift. On HST, the infrared channel of Wide Field Camera 3 will offer a big 
advance for faint galaxy surveys, offering a field of view 5.5 x larger than that 
of NICMOS, with 1.7x better pixel sampling. Deep 1.0-1.6 /jm observations 
of HDF-sized regions will become routine, enabling wholesale studies of the 
rest-frame optical morphologies, luminosities and colors of galaxies at z < 3, 
and color-selected surveys for galaxies at z > 6. 

Despite the promises usually made for NIR surveys (see §1), our census 
of the stellar content of galaxies at z > 2 is still fundamentally limited by the 
wavelengths at which we can observe. The H and iiT-bands measure rest- 
frame B and y-band light at z = 3, leaving large uncertainties on estimates 
of stellar mass, ages, and other such parameters, while brave first attempts at 
3-7 fj,m do not go deep enough to detect galaxies at z > 2 |3^||. The SIRTF 
Infrared Array Camera (IRAC), observing at 3.6-8.0 /xm, can measure rest- 
frame X-band light from galaxies out to z w 3 and Ao > 1 /im emission out 
to z = 7. Extremely deep exposures will be required, however, to detect ordi- 
nary galaxies at such large redshifts. We will be carrying out a SIRTF Legacy 
Program, the Great Observatories Origins Deep Survey (GOODS), which will 
push observations at 3.6-24 /im to their limits in two fields (the HDF-N and 
Chandra Deep Field South) totaling approximately 330 arcmin 2 . The survey 
goal is to provide multiwavelength data suitable for tracing the mass assem- 
bly history of galaxies and their energetic output from star formation and 
AGN activity out to the highest accessible redshifts. The SIRTF data, along 
with extensive supporting observations from ESO and other facilities, will 
be distributed to the community, providing a rich archive for research and a 
pathfinder to future work with NGST. 

MD would like to thank the conference organizers for hosting this im- 
portant and timely meeting. This work was supported by NASA grant GO- 
07817.01-96A. 
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